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1. INTRODUCTION 

Electromyography (EMG) signal contains rich muscle information that can be used in clinical and 
rehabilitation application. The potential of EMG signal in myoelectric control has been widespread since last 
two decades [1]. EMG signal recorded from a contracting muscle not only measures the time detection of 
muscle activation but also provides electrical signs of muscular behavior [2]. Recently, the analysis of EMG 
signal using a powerful signal processing technique has become the attention of the researchers. 

In biomedical signal processing, short time Fourier Transform (STFT), wavelet transform (WT) and 
empirical decomposition mode (EMD) are frequently used [3]-[5]. In the previous research, it has been found 
that WT outperformed other time-frequency methods in discriminating EMG patterns [3],[6]. WT exhibits 
good time resolution at high frequency and good frequency resolution at low frequency components [7]. 
In general, WT can be categorized into discrete and continuous form. In continuous wavelet transform 
(CWT), the wavelet transformation changes continuously. On one side, discrete wavelet transform (DWT) 
decomposes the signal into multiresolution coefficients using high pass and low pass filters. 
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Most studies to date indicated the performances of CWT and DWT were depending on the selection 
of a mother wavelet function [3],[8]-[10]. In the past studies, Kakoty et al. [8] investigated the best mother 
wavelet in DWT and CWT at different scale and decomposition level. The authors recommended the 
Gaussian and Symlet 8 to be the optimal mother wavelets in CWT and DWT, respectively. Phinyomark et al. 
[11] suggested that the use of DWT with the Daubechies 7 and 8 to ensure higher classification accuracy. 
Omari et al. [6] studied four mother wavelet functions at four different decomposition levels. The authors 
reported Symlet 4 offered the low classification error rate. Previous studies showed that the analysis of best 
mother wavelet in WT is critically important, leading to the optimum classification performance. However, 
the selection of mother wavelet is remains challenging in many areas. 

The best mother wavelet is mostly subject independent, which means different mother wavelet 
offers different kind of performance on different subject. In addition, previous works mostly focus on four to 
eight mother wavelets in the classification of EMG signals, which is insufficient. Moreover, the performance 
of mother wavelet at different scale and decomposition level provide significant difference in classification 
performance. It is obvious that the analysis of the mother wavelet in CWT and DWT is remain insufficient 
and unclear in EMG pattern recognition. Therefore, this study aims to evaluate the best mother wavelet in 
CWT and DWT by employing a large number of mother wavelet functions with different scale and 
decomposition level, respectively. 

This paper presents a detail study of the selection of mother wavelet in DWT and CWT. 14 mother 
wavelets of DWT and 12 mother wavelets of CWT at three different decomposition levels and scales are 
investigated, respectively. Two popular features mean absolute value (MAV) and wavelength (WL) are 
extracted from each wavelet coefficient for performance evaluation. The multiclass support vector machine 
(SVM) is used to classify EMG signal since it offers better performance in previous work [8],[12]. Finally, 
the best mother wavelet of DWT and CWT that offer the best classification performance are pointed. 


2. MATERIAL AND RESEARCH METHOD 
2.1. EMG data collection 

This study was performed on ten healthy subjects (8 males and 2 females) with mean age of 28.6 
(a=9.7) years. Each subject provided informed consent to participate in the experiment. Additionally, all 
subjects were free from neurological and muscular disorder. Two wearable EMG devices named Shimmer 
(Shimmer3 Consensys EMG Development Kits) with standard setting were used in data collection. The 
resolution was set at 24 bits with a gain of 12. The EMG signal was gathered from four useful hand muscles 
namely extensor digitorum (ch1), flexor carpi radialis (ch2), extensor carpi radialis longus (ch3) and flexor 
carpi ulnaris (ch4) with two reference electrodes at the elbow. The signal was sampled at 1024 Hz and band- 
pass filtered between 20 and 500 Hz. The skin was shaved and cleaned with alcohol pad before the electrode 
placement. The surface electrodes with 30 mm diameter were used and the inter-electrode distance was set at 
20 mm to reduce the crosstalk. The bipolar electrode configuration was shown in Figure 1. 


Figure 1. Electrodes configuration 


Subject was seated comfortably on a chair with the hand in neutral position. The surface EMG 
signals were recorded as the subject performed ten different hand movements including thumb flexion (M1), 
thumb extension (M2), wrist flexion (M3), wrist extension (M4), making a fist (M5), pinch index to the 
thumb (M6), pinch middle to the thumb (M7), pinch ring to the thumb (M8), pinch little to the thumb (M9) 
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and rest (M10). The experiments consisted of ten trials. Within each trial, the subject was asked to perform 
ten different hand movements for 5 s each, followed by a resting state of 4 s. Moreover, a resting period of 1 
min was introduced at the end of trial to avoid mental and muscle fatigue. The resting state was removed 
before data segmentation. 

A recent report of real time EMG application indicated that the optimal window length was ranging 
from 150 to 250 ms to balance the controller delay and classification error rate [13]. Additionally, an 
overlapped windowing technique was introduced to produce better classification accuracy in EMG pattern 
recognition [14]. In this work, the EMG data were divided into 250 ms window (256 samples) with 50% (128 
samples) overlapped. In total, a data matrix of 39 segments x 256 samples x 4 channels were obtained from 
each movement from each subject. 

Figure 2 shows the flow diagram of the proposed recognition system. In the first stage, the raw 
EMG signals are collected and segmented. Next, MAV and WL features are extracted from CWT and 
reconstructed DWT coefficients at different scale and decomposition level using different mother wavelet, 
respectively. In the final stage, the SVM is used to recognize the EMG signals of ten different hand 
movements. 


Raw EMG data 


Twelve mother 
wavelet 


Fourteen mother 
wavelet 


Scale, 2° 
(8, 16 and 32) 


Decomposition level 
(2, 4 and 6) 


Reconstruction 


Feature Extraction (MAV, WL) 


Classification (SVM) 


Hand Movement Types 


Figure 2. The flow diagram of the proposed recognition system 


2.2. Wavelet Transform 

Wavelet transform (WT) is a powerful mathematical tool that is successful in the analysis of bio- 
signal including EMG signal. WT offers high frequency resolution for low frequency component and good 
time resolution for the high frequency component [13]. Generally, WT can be categorized into continuous 
and discrete forms. Continuous wavelet transform (CWT) decomposes the signal based on the dilations and 
translations of a single mother wavelet function. CWT is more consistent and efficient because it provides 
localization time-frequency information without down-sampling [11]. Additionally, CWT is continuous in 
term of shifting and it gives useful time-frequency information [15]. CWT can be defined as: 


CWT. (s,b) = J xOy, (t)dt (1) 


where x(t) is the input signal and y;,(t) is the transformation of the mother wavelet function. 
The transformation can be expressed as: 


A Detail Study of Wavelet Families for EMG Pattern Recognition (Jingwei Too) 


4224 A ISSN: 2088-8708 


y o=- =) 
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where s is the scaling parameter, b is referred to the translation parameter and Yt) is called mother wavelet. 
The variables s and b provide the time scaling and shifting operation, respectively [16]. By using equation 1 
and 2, CWT can be computed as: 


(2) 


CWT(s,b) = 7 j ov ( Ja 
og : (3) 


Figure 3 demonstrates the scalogram of CWT at scale 32 using Mexican hat wavelet. The yellow 
areas represent higher amplitude at each scale. In turn, dark blue areas refer to low amplitude. 
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Figure 3. Scalogram of continuous wavelet transform at scale 32 using Mexican hat 


Discrete wavelet transform (DWT) is derived from CWT [17]. DWT is more widely used because it 
offers low computation cost [11]. In DWT, the signal is decomposed into the approximation and detail 
coefficient which involves the change of sampling rate [18]. The decomposition of DWT comprises of two 
digital filters, which are high-pass and low-pass filters. The low-pass and high-pass filter down-sample the 
input signal and provide the approximation, A and detail, D, respectively [11],[19]. For each decomposition 
level, the filters down-sample the signal by the factor of 2. The first level of decomposition is defined as: 


D[n] = $ xtk]-A[2n—k] 
$ (4) 


A{n] =) xik]: g[2n—k] 
7 (5) 


where x/k] is the input signal, D/n] is referred to the detail, D; and Afn] is the approximation, A;. 
The decomposition process is repeated until the desired final level is achieved. In the previous research, each 
coefficient subset was reconstructed to obtain more reliable EMG signal part, resulting in better classification 
accuracy [3],[13]. Therefore, the inverse wavelet transform is used to reconstruct each wavelet coefficient 
into more effective subset, namely, estimated approximation, rA and estimated detail, rD. For example, the 
estimated subset rD3 is obtained by performing the inverse wavelet transform on third-level detail, D3. 
The wavelet reconstruction of estimated detail (rD;-rD6) and estimated approximation (rA;-rA6) were shown 
in Figure 4. 
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Figure 4. Wavelet reconstruction of DWT at sixth decomposition level using Symlet 4 


2.3. Mother Wavelet Selection and Evaluation 
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Recent studies indicated WT has been recognized as one of the best time-frequency method in 
biomedical signal processing [3],[18],[20]. However, the performance of WT is mostly depending on the 
mother wavelet function. The selection of mother wavelet is remained challenging in many areas. Therefore, 
this work aims to evaluate the best mother wavelet in DWT and CWT for EMG signal processing. In this 
study, 14 mother wavelets in DWT and 12 mother wavelets in CWT are investigated. Table 1 is a lookup 
table of the mother wavelet used in CWT and DWT. It is worth noting different scale and decomposition 
level in CWT and DWT provide different property. For this reason, the performance of the mother wavelet at 
the scale 8, 16, 32 and decomposition level of 2, 4 and 6 are examined. 


Table 1. Mother wavelet of CWT and DWT used in this study 


CWT DWT 
1 Haar Haar 
2 Daubechies 2 Daubechies 2 
3 Daubechies 4 Daubechies 4 
4 Daubechies 6 Daubechies 6 
5 Symlet 2 Daubechies 8 
6 Symlet 4 Daubechies 10 
7 Symlet 6 Symlet 2 
8 Morlet Symlet 4 
9 Mayer Symlet 6 
10 Mexicanhat Symlet 8 
11 Gaussian 2 Coiflet 2 
12 Gaussian 4 Coiflet 3 
13 - Coiflet 4 
14 - Coiflet 5 


2.4. Feature Extraction using Wavelet Transform 


Feature extraction is an essential step to reduce the dimensionality and extract the useful information 
from the signal. In this work, wavelength (WL) and mean absolute value (MAV) are extracted from each 


wavelet coefficient. MAV and WL can be expressed as [6]: 
L 


1 
MAV =—), 


n=l 


X, 


n 


(6) 
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L-l 
WL = pare = Xa 
n=l 


(7) 


where X, is the input signal and L is the length of signal. 


2.5. Support Vector Machine 

Support vector machine (SVM) has been proved to be an outstanding supervised machine learning 
method in EMG pattern recognition [14]. Moreover, SVM has shown its superiority, especially for non-linear 
and high dimensional pattern recognition [21]. SVM maps the predictors onto a high dimensional space by 
using the concept of hyperplane partition for the data [22]. Some drawbacks of SVM are the complexity of 
the selection of kernel function and the longer computation time [14]. A previous study reported that radial 
basis function (RBF) was the best kernel function because it gave a higher classification performance [6]. In 
this regard, SVM with RBF kernel function is applied and it can be defined as: 


be 2 
K(x,x;) = oxo oat 


20° 


(8) 


where x-x; is the Euclidean distance between feature vectors and ais the kernel parameter. 


3. RESULTS AND ANALYSIS 

In this work, 10-fold cross validation is applied in the classification of EMG signals. The data is 
separated into 10 equal parts. Every part takes turn to test and the remaining parts are used in training phase. 
In the first part of the experiments, 14 mother wavelet functions in DWT at the three different decomposition 
level are evaluated. Table 2 outlines the mean classification accuracy of 14 mother wavelets of DWT at a 
decomposition level of 2, 4 and 6 across ten different subjects. From the results, the mean classification 
accuracy is found to be above 97% for all 14 mother wavelet functions in both WL and MAV feature sets. 
Additionally, MAV has shown to be an effective and reliable feature because it offers better performance in 
discriminating EMG patterns. By employing MAV feature, it is obvious that the highest classification 
accuracy is obtained by Symlet 4 (98.74%), followed by Daubechies 4 (98.72%) at the second decomposition 
level. On the other hand, Coiflet 3 outperforms other mother wavelets with a mean classification accuracy of 
98.49% at the fourth decomposition level when WL is used. From the analysis, Symlet 4 and Daubechies 4 at 
the second decomposition level are found to be the most suitable mother wavelet in DWT. 


Table 2. Classification Accuracy (mean + STD) of 14 Mother Wavelets of DWT at Three Different 


Decomposition Level Across Ten Subjects 
Classification performance (%) 


Classification performance (%) 


Mother wavelet Mother wavelet 


WL MAV WL MAV 
Level 2 97.90 + 1.02 98.43 + 0.88 Sym 4 Level 2 98.09 + 0.92 98.74 + 0.66 
Haar Level 4 98.00 + 0.90 98.28 + 0.80 Level 4 98.36 + 0.78 98.53 + 0.67 
Level 6 97.28 + 0.94 97.64 + 0.85 Level 6 97.50 + 0.79 97.67 + 0.81 
Level 2 97.97 + 1.01 98.63 + 0.68 Sym 6 Level 2 98.18 + 0.87 98.67 + 0.76 
Db2 Level 4 98.31 + 0.72 98.44 + 0.70 Level 4 98.39 + 0.67 98.55 + 0.68 
Level 6 97.32 + 0.94 97.45 + 0.78 Level 6 97.58 + 0.84 97.65 + 0.87 
Level 2 98.08 + 0.88 98.72 + 0.67 Sym 8 Level 2 98.19 + 0.87 98.70 + 0.71 
Db4 Level 4 98.36 + 0.78 98.56 + 0.68 Level 4 98.45 + 0.69 98.57 + 0.72 
Level 6 97.36 + 0.91 97.55 + 0.82 Level 6 97.67 + 0.87 97.74 + 0.85 
Level 2 98.23 + 0.90 98.65 + 0.73 Coif 2 Level 2 98.10 + 0.90 98.69 + 0.70 
Db 6 Level 4 98.48 + 0.63 98.53 + 0.69 Level 4 98.34 + 0.79 98.52 + 0.66 
Level 6 97.48 + 0.88 97.52 + 0.85 Level 6 97.59 + 0.91 97.70 + 0.77 
Level 2 98.20 + 0.90 98.69 + 0.71 Coif 3 Level 2 98.18 + 0.90 98.69 + 0.71 
Db 8 Level 4 98.44 + 0.67 98.59 + 0.67 Level 4 98.49 + 0.70 98.62 + 0.62 
Level 6 97.57 + 0.74 97.60 + 0.82 Level 6 97.56 + 0.85 97.71 + 0.73 
Level 2 98.17 + 0.94 98.70 + 0.68 Coif 4 Level 2 98.22 + 0.93 98.70 + 0.70 
Db 10 Level 4 98.48 + 0.66 98.62 + 0.60 Level 4 98.42 + 0.74 98.56 + 0.68 
Level 6 97.43 + 0.91 97.49 + 0.88 Level 6 97.61 + 0.83 97.71 + 0.72 
Level 2 97.97 + 1.01 98.63 + 0.68 Coif 5 Level 2 98.23 + 0.88 98.70 + 0.71 
Sym 2 Level 4 98.31 + 0.72 98.44 + 0.70 Level 4 98.45 + 0.72 98.56 + 0.59 
Level 6 97.32 + 0.94 97.45 + 0.78 Level 6 97.50 + 0.86 97.59 + 0.85 
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In the second part of the experiments, 12 mother wavelets of CWT are studied. Table 3 
demonstrates the mean classification accuracy of 12 mother wavelets of CWT at scale 8, 16 and 32 for ten 
different subjects. At scale 8, Gaussian 2 and 4 exhibit the highest classification accuracy of 98.42% using 
WL and MAV feature sets, respectively. However, their performance did not show much improvement at a 
higher scale. At scale 16, it has been found that the Symlet 6 achieves the best classification accuracy of 
98.56%, followed by Symlet 4, 98.53% when MAV is used. For instance, the Mexican hat has shown its 
superiority at scale 32 with the best mean classification accuracy of 98.64% in WL feature set. Unfortunately, 
MAV shows the decrement in classification performance at scale 32. This shows that MAV feature set is not 
suitable for high scale wavelet function in CWT. As a result, the most suitable mother wavelet in CWT are 
Mexican hat at scale 32 and Symlet 6 at scale 16. 


Table 3. Classification Accuracy (mean + STD) of 12 Mother Wavelets of CWT at Three Different Scale 
Across Ten Subjects 


Classification performance (%) Classification performance (%) 


Mother wavelet Mother wavelet 


WL MAV WL MAV 
Scale 8 97.70 + 0.96 98.00 + 1.08 Sym 6 Scale 8 98.05 + 0.86 98.17 + 0.97 
Haar Scale 16 98.38 + 0.92 98.31 + 0.96 Scale 16 98.48 + 0.81 98.56 + 0.74 
Scale 32 98.51 + 0.76 98.19 + 0.79 Scale 32 98.49 + 0.72 98.35 + 0.73 
Scale 8 97.88 + 1.01 98.06 + 1.13 Morl Scale 8 98.00 + 0.83 98.07 + 0.83 
Db 2 Scale 16 98.44 + 0.88 98.42 + 0.86 Scale 16 98.40 + 0.86 98.40 + 0.82 
Scale 32 98.50 + 0.72 98.29 + 0.73 Scale 32 98.34 + 0.74 98.26 + 0.78 
Scale 8 97.99 + 0.92 98.13 + 1.03 Meyr Scale 8 98.06 + 0.95 98.13 + 0.95 
Db4 Scale 16 98.46 + 0.90 98.47 + 0.78 Scale 16 98.45 + 0.84 98.49 + 0.75 
Scale 32 98.45 + 0.76 98.27 + 0.76 Scale 32 98.36 + 0.79 98.29 + 0.81 
Scale 8 97.94 + 1.01 98.08 + 1.08 Mexh Scale 8 98.36 + 0.82 98.15 + 0.79 
Db6 Scale 16 98.41 + 0.93 98.46 + 0.78 Scale 16 98.08 + 0.76 97.49 + 0.81 
Scale 32 98.36 + 0.75 98.27 + 0.76 Scale 32 98.64 + 0.66 96.26 + 1.00 
Scale 8 97.88 + 1.01 98.06 + 1.13 Gaus 2 Scale 8 98.42 + 0.83 98.35 + 0.84 
Sym 2 Scale 16 98.44 + 0.88 98.42 + 0.86 Scale 16 98.28 + 0.76 98.00 + 0.77 
Scale 32 98.50 + 0.72 98.29 + 0.73 Scale 32 98.50 + 0.67 97.01 + 0.87 
Scale 8 98.03 + 0.87 98.18 + 0.99 Gaus 4 Scale 8 98.39 + 0.93 98.42 + 0.83 
Sym 4 Scale 16 98.48 + 0.83 98.53 + 0.74 Scale 16 98.48 + 0.70 98.42 + 0.71 
Scale 32 98.52 + 0.69 98.34 + 0.74 Scale 32 98.31 + 0.69 97.80 + 0.77 


In the final part of the experiments, the paired two-tail t-test is used to measure the statistical 
difference between the classification performances of WL and MAV features when different mother wavelet 
function is used. Table 4 and 5 outline the result of t-test of the classification performance obtained from 
DWT and CWT across ten subjects. In t-test, the null hypothesis is rejected if the p-value is less than 0.05. 
This shows that there is a statistical difference between WL and MAV feature sets. 

From Table 4, the results of the WL and MAV are statistical difference for all wavelet functions at 
the second decomposition level. At fourth decomposition level, the p-value illustrates that the Daubechies 6 
and Coiflet 5 show no significant difference when WL versus MAV. At sixth decomposition level, only 
Haar, Daubechies 4 and Symlet 4 exhibit the significant difference. From Table 5, Haar, Symlet 4 and 
Mexican hat show significant difference in scale 8. Additionally, at scale 16, only Mexican hat, Gaussian 2 
and Gaussian 4 obtain p-value lower than 0.05. Moreover, other than Daubechies 6 and Symlet 6 exhibit 
significant differences between the classification performance of WL and MAV at scale 32. 


Table 4. Result of t-test of the Classification Performance between MAV and WL using DWT 


p — value 
Mothier Wavelet Level 2 Level 4 Level 6 
Haar 0.0007 0.0007 3E-05 
Db2 0.0012 0.0195 0.0521 
Db4 0.0006 0.0087 0.0031 
Db6 0.0070 0.3754 0.2340 
Db8 0.0037 0.0185 0.6085 
Db 10 0.0020 0.0036 0.3163 
Sym 2 0.0012 0.0195 0.0521 
Sym 4 0.0009 0.0057 0.0138 
Sym 6 0.0008 0.0046 0.0380 
Sym 8 0.0007 0.0289 0.1081 
Coif 2 0.0012 0.0178 0.0854 
Coif 3 0.0031 0.0109 0.0625 
Coif 4 0.0031 0.0157 0.0504 
Coif 5 0.0010 0.0860 0.0807 
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Table 5. Result of t-test of the Classification Performance between MAV and WL using CWT. 


p — value 
Motherwavelet Scale 8 Scale 16 Scale 32 
Haar 0.0377 0.2676 0.0003 
Db2 0.0539 0.8141 0.0003 
Db4 0.1104 0.8478 0.0104 
Db6 0.0525 0.3855 0.0670 
Sym 2 0.0539 0.8141 0.0003 
Sym 4 0.0409 0.5256 0.0037 
Sym 6 0.0625 0.2172 0.0526 
Morl 0.0635 0.9162 0.0050 
Meyr 0.1266 0.3865 0.0207 
Mexh 3E-05 7E-07 2E-06 
Gaus 2 0.2864 TE-05 TE-07 
Gaus 4 0.5683 0.0334 5E-05 


4. CONCLUSION 

In this study, the usefulness of the mother wavelet function in DWT and CWT has been 
investigated. Two popular features, WL and MAV are extracted from the wavelet coefficients as the input to 
the SVM. In CWT, the Mexican hat at scale 32 and Symlet 6 at scale 16 are suggested to be the optimal 
mother wavelet selection for the classification of EMG signals. On the other hand, the reconstructed DWT 
coefficient with Daubechies 4 and Symlet 4 at second decomposition level are recommended to be used in 
EMG pattern recognition. The experimental results indicated DWT not only offered low computation cost, 
but also yielded a high classification accuracy. As compared to CWT, DWT is more approariate to be used in 
rehabilitation and clinical application. 
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